Direct estimation of single- and two-qubit Hamiltonians and relaxation rates 



oo 

o 

O 

(N 

Oh- 
< 

(N 



43 ■ 

a,; 

H— > ■ 

c ■ 

ctf ■ 
=5 

cr 



> 

p 
O 

o 



X 



M. Mohseni, 1 A. T. Rezakhani, 2,3 and A. Aspuru-Guzik 1 

'Department of Chemistry and Chemical Biology, Harvard University, 12 Oxford St., Cambridge, MA 02138, USA 
2 Center for Quantum Information Science and Technology, and Departments of Physics and Chemistry, 
University of Southern California, Los Angeles, CA 90089, USA 
3 Institute for Quantum Information Science, University of Calgary, Alberta T2N 1N4, Canada 

We provide a novel approach for characterization of quantum Hamiltonian systems via utilizing a single mea- 
surement device. Specifically, we demonstrate how external quantum correlations can be used for Hamiltonian 
identification tasks. We explicitly introduce experimental procedures for direct estimation of single- and two- 
qubit Hamiltonian parameters, and also for simultaneous estimation of transverse and longitudinal relaxation 
rates, using a single Bell-state analyzer. An advantage of our method over the earlier approaches is that it has a 
built-in feature which makes it suitable for partial characterization of Hamiltonian parameters. 

PACS numbers: 03.65.Wj, 03.67.Lx 



I. INTRODUCTION 

Characterization of quantum dynamical systems is one of 
the most fundamental problems in quantum physics, and lies 
at the heart of quantum information processing [jj] and coher- 
ent control iH]. There are a few known methods to achieve 
this goal, such as standard quantum process tomography 0], 
ancilla-assisted process tomography J4J], and direct character- 
ization of quantum dynamics (DCQD) JH 0, 0]. Since the 
required physical resources grow exponentially with the num- 
ber of degrees of freedom, all quantum process tomography 
schemes are in principle inefficient. However, in various phys- 
ical situations a full characterization of the quantum dynam- 
ical superoperator is not always necessary, as sometimes the 
information about relevant physical quantities could be related 
to only a polynomial number of parameters in the system size 
id. This is indeed the case when: (i) important physical prop- 
erties of a quantum system can be directly associated only to a 
subset of certain superoperator elements, (ii) a priori knowl- 
edge exists about general properties of quantum dynamics, 
and (iii) neglecting some elements will lead to small system 
characterization errors. 

The task of Hamiltonian identification, as a characteri- 
zation of quantum systems, is of paramount importance in 
quantum physics and chemistry. In particular, it is required 
for monitoring or controlling performance of noisy single- 
and two-qubit quantum gates/devices in quantum information 
processing Ql |9[]. For various physical systems, a generic 
form of the Hamiltonian can be guessed from general physi- 
cal/engineering considerations or observations. However, one 
still needs to estimate the Hamiltonian parameters for a given 
quantum system in order to study the internal dynamics of the 
system and also to investigate the exact form and the strength 
of a potential system-bath coupling. 

Identification of time-independent (or piece-wise constant) 
Hamiltonians along with the estimation of errors have al- 
ready been studied for both single-qubit and two-qubit cases 
iflol 1 1 ill . Characterization of a single-qubit Hamiltonian is 
achieved via determination of the measurement results using 
a single fixed readout process which is a periodic function of 
time. Through Fourier analysis of this signal and other re- 



lated techniques, identification is reduced to finding the (rel- 
ative) location of peaks and their heights of the Fourier spec- 
trum fioh . Similarly, in the two-qubit case, Hamiltonian pa- 
rameters are obtained through some entanglement measure- 
ments sampled in many time points, and then a Fourier anal- 
ysis determines the parameters fTlll . There are also more 
general Hamiltonian identification schemes which employ 
closed-loop learning control techniques along with efficient 
and improved global laboratory data inversion for identifica- 
tion 111211 . These techniques are useful particularly when one 
has access to tailored control fields (e.g., shaped laser pulses) 
while the measurements are being performed. A fundamental 
relevant question is how one can exploit external quantum cor- 
relations in order to enhance identification of quantum Hamil- 
tonian systems. This is the subject we address in this work. 

In this paper, we introduce an analytical method for direct 
characterization of important classes of Hamiltonians. This 
method is based on a newly proposed DCQD scheme @]. In 
particular, we demonstrate how to estimate all parameters of a 
general time-independent single-qubit Hamiltonian and two- 
qubit (isotropic or anisotropic) exchange Hamiltonian. A dis- 
tinctive feature of our method is that, when only some partial 
knowledge about the system is requested, it does not require 
Fourier analysis of the experimental data. In principle, this 
obviates the need for long sampling times and in turn offers 
more controllability for the related estimation process. Our 
Hamiltonian identification method is applicable to quantum 
systems enabling two-body measurements, due to the fact that 
DCQD requires Bell-state measurements (BSMs). The re- 
quired BSM can be in principle achieved in linear optical sys- 
tems via postselections lfl3ll or h ype rentanglement (3, and 
also in trapped ions lfl5ll (see Ref. lllotl for a deterministic, pro- 
grammed generation of "ultralong lifetime" Bell-states), and 
optical lattices lfl7ll . In solid-state systems, several schemes 
for controllable two-body interactions have been proposed 
d [H d HJ IH HI |2J, with the state of the art exper- 
imental realization in semiconductor quantum dots l24ll and 
superconducting flux qubits l25ll . 

The evolution of a rf-dimensional quantum system (open 
or closed) with initial state p can, under some natural as- 
sumptions, be expressed in terms of a completely positive 
quantum dynamical map £ t , which can be represented as: 
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£ t{p) = YjijXij{t)viPv)- The positive matrix x = [Xmn] 
encompasses all information about the dynamics, relative to 
the fixed operator basis set {er,„}, where tr(trj n cr n ) = d5 m n- 
The theory of DCQD determines elements of x matrix J5t] by 
relating them to measurement results more directly than the 
other existing schemes. The main idea of DCQD is based 
on quantum error detection theory in which by preparation of 
suitable states and measurement of their stabilizers and nor- 
malizers partial information about errors can be obtained. The 
required stabilizer and normalizer measurement can be phys- 
ically realized with a single Bell-state analyzer. Table U sum- 
marizes the scheme for the single-qubit case. 



H. IDENTIFICATION OF SINGLE-QUBIT 
HAMILTONIANS 

Let us consider the cases that quantum dynamics is 
generated by a time-independent Hamiltonian, £t{p) = 
U{tYp U(t), where U(t) = e~ ltH (h = 1), one obtains 
Xmn = a m a n , where U{t) = ^ m fl m (i)ff ra and H is 
the Hamiltonian of the system. Since an energy shift is al- 
ways possible, we only consider traceless Hamiltonians. In 
the single-qubit case, where H = J ■ <r, with the choice 
of {<r = (a x , a y , cr z ), 1} as the operator basis, we have: 
Xoo = c 2 , Xaa = s 2 J 2 , xoa = iscJ a , and Xa/3 = s 2 J a jp, 
in which J = J J (J — \\J\\), c = cos(Ji), s = sin(Jt) and 
a, 13 = x,y, z. 

In order to find the real vector J, according to the DCQD 
theory, we can choose different experimental configurations 
(measurement settings) depending on our a priori informa- 
tion about the Hamiltonian. If the signs of the Hamiltonian 
parameters (i.e., the components of J) are already known, we 
can determine J and the absolute values of the components, 
| J a \, in a single experimental configuration. First, we prepare 
a maximally entangled state between the qubit of interest A, 
and an ancilla, B, as \$ + )ab = ([00) + \11))/V2. Next, the 
system evolves under H for a duration of time t. Finally, we 
perform a BSM represented by the four projection operators 
P$± and P^±. The probabilities of obtaining these outcomes 
are found as tr[Pj£ t (p)] = c 2 and tr[P a £ t (p)] = s 2 J 2 , where 
the projection operators Pi and P a (for a = x,y,z) corre- 
spond to the states $ + , and $ _ , respectively. Thus, 



TABLE I: Input states and measurements for direct characterization 
of single-qubit X- Here |$+) = a|00) + P\U), \$t}x( y ) = a\ + 
+>*(») +P\~ -)x(y) where \a\ ^ # and Im(a/3) + 0, and 
{|0), |1)}, {|±)a:}, {|±) H } are eigenstates of the Pauli operators a z , 
a x , and a y . is the projector on the Bell state |$ + ), and similarly 
for the other projectors. 



input state 


Bell-state measurement 


OUtpUt Xmn 


l*+) 




X00, Xll, X22, X33 


i*i> 


p s+ ±p^-,p q ,+ ±p <[ ,- 


X03, X12 




p s+ ±p^+,p i> - ±py- 


X01 i X23 




p s+ ±p^_,p # _ ±P (P+ 


X02, X13 



we have 

J a 2 - tT[P a £ t (p)]/(l - tr[Pr£ t (p)]). (1) 

The diagonal elements of the superoperator give the absolute 
values of the unknown parameters J a . Equation (|TJ bears this 
interesting result that measurements at a single time-point t 
are in principle (ignoring the inherent issue of statistical er- 
rors) enough to obtain (| J x \, \ J y \, \ J z \), When the relative 
signs are already known, this uniquely identifies the reference 
frame of the Hamiltonian. 

In order to obtain J, we are required to estimate the fre- 
quency of the function cos 2 ( Jt) = tr [Pj£ t (p)]. The theory of 
signal processing and discrete Fourier analysis state that one 
generally needs to perform many time samplings to obtain fre- 
quencies. By the Nyquist criterion, the sampling frequency 
fs = 1 /ts must be bounded below by half of the frequency 
of the original signal, i.e., f,s > J, to reduce the inherent 
aliasing |2fl]. In Refs. |10UllH one can find more detailed 
analysis of these issues and how to read J from experimen- 
tal data. Specifically, in Ref. ITlVb)! an interesting method of 
ensemble measurements in sample points has been introduced 
that can reduce the statistical error in inference. 

In the more general case, to fully characterize the real 
vector J we need to consider a different strategy and per- 
form two measurements for the off-diagonal elements of X- 
According to DCQD, these two experimental configurations 
are sufficient to determine the diagonal of the superopera- 
tor, xu f° r * £ {0,1,2,3}, and four off-diagonal param- 
eters Im(xoi)' and Re(xjk), for any two sets of values of 
{i,j,k}, {i',j',k'} e {1,2,3} such that i ^ i', j ^ k ^ i 
and j' k' ^ i' . For example, by preparation of a non- 
maximally entangled state |<&+) = a|00) + /3|11) (Table[I]i 
and performing a standard BSM, we can obtain the following 
equations: 

X00 +X33 =P+, Xll+X22=P-, 

a(xoo - X33) + Mrn(xo3) = c + , 
a(xn ~ X22) - 6Re(xi2) = c_, 

withp± = tT[P ±1 £ t (p)}, a = tr(A/», b = 2itr(a^Afp), and 
c± = P±tr(A/p±i), where p = |$+)(*+|, N = a^trf, 
P + i = P 9 + + P $ -, P_! = P^+ + P 9 - and p ±1 = 
P±i£t(p) P±i / te[P±i£t(p)]- In the other experimental con- 
figuration, we prepare a nonmaximally entangled state 
and perform another standard BSM to obtain a similar set of 
equations in the {Izb)^} basis (here and also for {|±) y } ba- 
sis, Af = erf erf ). Using these linearly independent equations 
we can determine diagonal elements of the superoperator, xa 
(i = 0,1,2,3) and four off-diagonal parameters Im(xo3), 
Im(xoi), R e (xi2), and Re(x23)- As we have shown above, 
the diagonal elements can be used to determine J and the ab- 
solute values \J a \. The relative signs of J x , J y , and J z can 
be found from the off-diagonal parameters above; so, we can 
identify J up to a global sign. This global sign is usually 
evident from the physical/engineered setup under considera- 
tion, e.g., from the direction of a global magnetic field for spin 
systems. In physical situations where this global sign cannot 
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be deduced from general physical considerations, we need to 
perform a third measurement that corresponds to character- 
izing Im(xo2) and Re(x3i) which completes our knowledge 
about an arbitrary (time-independent) single-qubit Hamilto- 
nian. The whole analysis is also applicable to the case of 
piece-wise constant Hamiltonians. In the following we dis- 
cuss another important example of single-qubit dynamics, al- 
though non-Hamiltonian, that shows how the DCQD estima- 
tion may provide advantage in estimation of dynamical pa- 
rameters in the Markovian regime. 



III. SIMULTANEOUS DETERMINATION OF Ti AND T 2 

Let us consider the so-called quantum homogenization pro- 
cess acting on a single-qubit density matrix p(0) for time 
t, where poo(0) = a and poi(O) = b in the {|0), |1)} ba- 
sis. This leads to the final state p{t) with poo(t) = (a — 
doo) exp(— i/Ti) + and poi(t) = bexp(-t/T 2 ), where 
doo characterizes the population of thermal equilibrium state, 
and the time-scales Ti and T2 (T2 ^ 2Ti) are longitudi- 
nal and transverse relaxation time-scales of the system, re- 
spectively 01 |27l]. The explicit form of xu elements are 
as follows: X oo(33) = [exp(-i/Ti) ± 2exp(t/T 2 ) + l]/4, 
X11 = X22 = [exp(t/T x ) + l]/4. 

Now we demonstrate that both T\ and T2 can always be 
estimated in a single ensemble measurement by using the 
DCQD scheme for estimating diagonal elements of \- We 
first prepare a Bell-state \& + )ab, and then let the qubit A in- 
teract with a thermalizing environment for a given time t. The 
outcomes of a BSM yield the following relations for T\ and 
T 2 : 

1/Ti = - In (2tr[P*+£ t (p)] + 2tr[P*-£ t (p)] - 1) ft, 
1/T a = -ln(tr[P* + £t(p)] -tr[P*-£(p)])/t. 

Ideally, these equations imply adequacy of single time-point 
measurements. That is, unlike the case of reading J, where 
time sampling is necessary and aliasing is inevitable, T\ and 
T2 can in principle be obtained through single time-point mea- 
surements. This feature could be utilized to reveal the non- 
Markovian nature of system-bath interaction; e.g., by observ- 
ing time variations in the estimated relaxation and dephasing 
rates beyond the the natural deviation due to statistical errors. 
Moreover, due to orthogonality of BSM outcomes, it is easy 
to unambiguously distinguish T\ from T2, unlike the approach 
presented in Ref. l3~lll . Traditionally, in order to measure 
the longitudinal and transverse relaxation times, one needs to 
measure two non-commutative observables (e.g., Pauli oper- 
ators a z and a x ) on two subensembles of identical systems, 
corresponding to magnetization vectors M z and M my parallel 
and perpendicular to a global magnetic field Bq. The num- 
ber of repetitions in each measurement is determined by the 
desired accuracy in the time sampling estimation of the relax- 
ation times associated with magnetizations M z and M xy l32tl . 



IV. IDENTIFICATION OF TWO-QUBIT EXCHANGE 
HAMILTONIANS 

In solid-state systems, it is often the case that each pair 
of qubits (AB) interact directly or effectively through an ex- 
change Hamiltonian H ex = ^2 a J a cr^a^, where J a s are 
the couplings of the two-qubit interaction (see also Ref. lfl7ll 
for the exchange interaction between neutral atoms in opti- 
cal lattices). The case of isotropic or Heisenberg interaction 
( J x = J y = J z ) is intrinsic to spin-coupled quantum dots, 
and donor atom nuclear/electron spins 1180 . This interac- 
tion is also important in the context of universal fault-tolerant 
quantum computing 111911 . The XY Hamiltonian (J x = J y , 
J z = 0) is the available interaction in quantum Hall sys- 
tems J20I1 . quantum dots/atoms in cavities I2H1 . and exciton- 
coupled quantum dots l22tl . The XXZ (J x = J y ^ J z ^ 0) 
interaction appears in the electrons in liquid-Helium quantum 
computing proposals i23ll . 

In the case of XYZ Hamiltonians, the x matrix has only 10 
nonzero independent elements Xmn, for m,n = 0,5, 10, 15. 
Similar to the case of the general single-qubit Hamiltonian, 
these diagonal elements contain information only about the 
absolute values | J a \s. In order to obtain information about 
the signs of J a s, we need to measure off-diagonal elements 
as well. However, in most physical/practical cases the signs 
of the terms in an exchange Hamiltonian are already known 
from some general properties of the system. For example, for 
many materials it is known whether below the phase transi- 
tion point they become ferromagnetic or anti-ferromagnetic — 
alternatively this information can be obtained for a given ma- 
terial by measuring its linear response to an applied magnetic 
field. In these cases, the Hamiltonian can be completely de- 
termined with a single ensemble measurement corresponding 
to the diagonal elements of the superoperator. 

Let us consider the important classes of isotropic and 
anisotropic exchange interactions. For these Hamiltonians the 
sign of J is known from the ferromagnetic property of the sys- 
tem. In fact, by definition J = E$ — Et (where E$ and Et 
are the energy of singlet and triplet states), is always negative 
for ferromagnetic materials. For example, for a two-electron 
system, the singlet state is the ground state of the system if 
J < 0. On the contrary, for anti-ferromagnetic materials, J is 
always positive which indicates that in the ground state spins 
tend to arrange themselves in the same direction. 

In order to determine the diagonal elements xu f° r Heisen- 
berg interaction between two electrons A\ and Ai, one can 
prepare a tensor product of maximally entangled states be- 
tween each electron and a pair of ancilla electrons (B\ and 
B2) such as \& + )a 1 b 1 \® + )a 2 b 2 - Then, the unknown Hamil- 
tonian H for the duration of t is applied, and a tensor 
product of BSMs acting on each pair AiBi is performed, 
where this operation can be represented by a tensor prod- 
uct of P $ +, Pj+, Pp-, Pj,- for i = 1,2, The joint prob- 
ability distributions of the BSMs are related to J through 
tr[Pj,+ P $ +£ t (p)] = c 6 + s e and tr[P^+ £ t {p)] = 

tr[P* r P* a -^(p)] = tT[P 9 -P^-S t {p)] = s 2 c 2 .~ Therefore, 
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we have: 

sin(2| J\t) = 2^tv[P^-P^-£ t (p)], (3) 

and similar relations hold for P^,+ P^+ and P^-P^- as well. 

In the case of anisotropic exchange interactions one can 
perform a similar Bell-state preparation and BSM as in the 
case of isotropic exchange, to obtain: 

sm(2\J x \t) = 2 y /tv[Py±P iS! ±£ t (p)], (4) 

cos(2| J z \t) = yf{ix[P 9t P 9t Stip)] - - 4)- (5) 

To read | J a |, one needs to have time samplings (i.e., ensemble 
measurements for many time-points) and follow the Fourier 
analysis based method sketched earlier. Therefore, having a 
priori knowledge about the ferromagnetic property of the sys- 
tem, one can identify the underlying exchange Hamiltonian. 

Note that the energy spectrum of H ex can be simply calcu- 
lated using the above relations and knowing the fact that Bell- 
states are the eigenkets of the exchange Hamiltonian. Eigen- 
values of H ex can be written as E = ± | J a | ± | Jp — J 1 | , where 
a 7^ (3 7^ 7 e {x, y, z}. We have already shown how to es- 
timate | J a \ for a — x,y, z. In order to find relative signs of 
any two other components, such as \ J y — J z \, the DCQD al- 
gorithm can be utilized by performing a single ensemble mea- 
surement that corresponds to measuring the off-diagonal ele- 
ment xo,5- For full characterization of an exchange Hamilto- 
nian without having any a priori knowledge about the signs of 
the coupling constants, one needs to measure the off-diagonal 
element xo.io too. Therefore, with a total of three ensem- 
ble measurements, corresponding to Xi,i> Xo,5> an d Xo,io, full 
characterization of H sx can be achieved. 



V. REMARKS ON PRECISION 

In a realistic estimation process, due to decoherence, lim- 
ited measurement or preparation accuracies (because of spe- 
cific device architecture or finite ensemble size), and other 
imperfections, some errors may occur (a generalization of 
the DCQD theory that addresses faluty preparations and 
measurements is underway and will be reported elsewhere 
l28tl ). These factors might affect the amount of physical 
resources required for a given accuracy of the estimation, 
hence some appealing features like sufficiency of single-time 
point measurements might be lost. For the case of ideal 
preparations/measurmens scenario in which single time-point 
measurements are in principle sufficient, errors scale up as 
1/VNe, where Ne is the number of repeated measurements. 
In the cases we need to perform time samplings, the error in 
the estimation of frequencies (and thus, Hamiltonian param- 
eters) is governed by the Nyquist criterion and the quantum 
shot-noise limit |2^1 . Let us consider iVg samples, for each of 
which we perform Ne measurements. Therefore, according 



to the quantum shot-noise limit, A/ ~ 1/(Nsts^/Ne), and 
the Nyquist criterion, fs_= l/ r S //2, we get: A/// ~ 
1/(N s VNe) (see Ref. 1111(b)]). That is, the average error in 
estimation of Hamiltonian parameters scales as 1/(Ns^/Ne~). 
In other words, for an error e, or with the number of digits of 
precision log(l/e), we need poly(l/e) more steps, which is 
common among all Fourier analysis based methods )T30tl - In 
this respect, our method does not provide an advantage over 
the one in Refs. ITHl — both methods provide similar accuracy 
scaleup. However, the advantage of our method lies in the 
built-in ability of the method for partial characterization. That 
is, there is a level of independency in the way different sets 
of parameters are related to measurements results. For more 
general discussions on partial characterization by DCQD see 
Refs. JH and for a very recent proof-of -principle experiment 
on this issue, using polarization and spatial degrees of free- 
dom of a single photon, see Ref. ^ . Moreover, in our method 
some of the parameters are related more directly to the mea- 
surements data, hence obviating the need to (a full) inversion 
in the first place. E.g., in the single-qubit case, we obtain 
|</ic,y,z| just by a very simple algebraic manipulation of the 
data. This feature is not necessarily available in the methods 
of Refs. fllll . because a Fourier analysis would be necessary 
even for extracting a partial information about the Hamilto- 
nian. 



VI. CONCLUSION AND DISCUSSION 

We have presented a theoretical approach for utilizing aux- 
iliary quantum correlations to perform Hamiltonian identifi- 
cation. In this method one can directly obtain full information 
about unknown parameters of time-independent of single- and 
two-qubit Hamiltonians without full quantum process tomog- 
raphy. In addition, we demonstrate that for a single qubit un- 
dergoing a generic Markovian homogenizing quantum map, 
both related relaxation times can be estimated simultaneously 
by utilizing a single Bell-state measurement. Furthermore, we 
illustrate how our prior knowledge about Hamiltonian systems 
can be exploited in order to reduce the required physical re- 
sources for identification tasks. In particular, we show that 
the required repeated measurements, associated to time sam- 
pling of data, can be reduced when we are interested in partial 
characterization of the Hamiltonian systems and also for es- 
timating relaxation rates. With the recent advent of various 
methods for generation of controllable entanglement, our pro- 
posed method may have near-term application to a variety of 
quantum systems/devices, such as in trapped ions, liquid-state 
NMR, optical lattices, and entangled pairs of photons. 
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(to M.M.), the Faculty of Arts and Sciences of Harvard Uni- 
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